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SUMMARY 


^  I  A  boundary -value-collocation  procedure  was  applied  in  conjunction  with  the 
Williams  stress  function  to  determine  values  of  the  stress -Intensity  factor  Kj 
for  single-edge  cracks  in  plate  specimens  subject  to  splitting  forces  applied 
close  to  the  crack  and  acting  transversely  to  it.  The  results  are  presented  in 
terms  of  the  dimensionless  quantity  Y  =  KjBH^/^/pa^  where  B  and  H  are  the 
specimen  thickness  and  half -depth,  a  is  the  effective  crack  length,  and  P  is 
the  applied  load.  The  results  are  practically  independent  of  the  ratio  of  ef¬ 
fective  specimen  length  to  specimen  half-depth  q/H,  when  this  ratio  is  not  less 
than  a/H  +  2,  and  are  then  in  excellent  agreement  with  those  derived  by  other 
investigators  from  compliance  measurements.  In  the  limit,  as  H/a  approaches 
zero,  the  value  of  Y  approaches  that  obtained  in  an  elementary  analysi^  ,v»5hich 
treats  the  specimen  as  a  pair  of  built-in  cantilever  beams .  |  / 1  ^  'J„./ 


INTRODUCTION 


Plate  specimens  of  the  general  configuration  shown  in  figirre  1  are  used  for 
plane  strain  fractirre  toughness  measurements  (ref.  l)  and  for  fatigue  crack 
propagation  rate  measirrements .  This  type  of  specimen  is  essentially  a  narrow 
plate  with  a  crack  extending  from  one  of  the  short  edges  along  the  longitudinal 
centerline.  The  forces  P  are  applied  close  to  the  crack  and  act  transversely 
to  it,  thus  tending  to  split  the  specimen  by  extending  the  crack.  Such  speci¬ 
mens  are  referred  to  in  this  report  as  single -edge -notch  specimens  subject  to 
splitting  forces.  They  are  sometimes  referred  to  elsewhere  as  crackline  loaded 
specimens,  or  as  double  cantilever  specimens.  The  length  of  the  crack  a  is 
taken  as  the  distance  from  the  end  of  the  crack  to  the  line  of  action  of  the 
load  P.  (All  symbols  used  in  this  report  are  defined  in  the  appendix.) 

A  first  approximation  to  the  strain  energy  release  rate  with  crack  exten¬ 
sion  for  such  a  specimen  is  readily  calculated  from  simple  beam  theory  by  re¬ 
garding  the  specimen  as  a  pair  of  built-in  cantilever  beams.  This  approach  was 
used  by  Gilman  in  a  study  of  the  mechanics  of  cleavage  (ref.  2).  In  terms  of 
the  stress-intensity  factor  Kj,  the  result  obtained  is  that  Kj  is  equal  to 


2  -yfz  Pa/H^/^  for  small  values  of  H/a,  where  H  is  the  depth  of  each  canti¬ 
lever  arm  (ref.  3).  A  better  approximation  was  obtained  by  Kipling,  Mostovoy, 
and  Patrick  (ref.  4),  who  modified  the  cantilever  beam  treatment  and  used  ex¬ 
perimental  compliance  measurements  to  establish  the  best  value  for  a  correction 
term. 


Because  of  the  potential  usefulness  of  this  type  of  specimen,  it  was  con¬ 
sidered  desirable  to  conduct  an  Independent  study  of  the  configuration  by  the 
boundary-collocation  procedure  used  in  earlier  studies  of  other  types  of  single¬ 
edge-notch  specimens  by  Gross,  Srawley,  and  Brown  (refs.  5,  6,  and  7).  It  was 
decided  that  this  study  should  also  include  an  evaluation  of  the  effect  of  a 
variable  not  considered  in  reference  4,  namely,  the  ratio  of  the  effective  spec¬ 
imen  length  to  the  effective  crack  length  q/a.  It  is  apparent  that  the  stress- 
intensity  factor  will  be  independent  of  q/a  for  sufficiently  large  values  of 
that  ratio,  but  it  is  of  practical  interest  to  know  how  large  a  ratio  is  suffi¬ 
cient. 

In  general,  the  stress -intensity  factor  Kj  will  depend  on  the  five  vari¬ 
ables  P,  a,  q,  H,  and  B  (fig.  l) .  It  is  convenient  to  express  the  relation 
between  Kj  and  four  of  these  variables  in  dimensionless  form  in  terms  of  the 

stress -intensity  coefficient  Y  =  KjBH^/^/pa,  which  is  a  function  of  the  ratios 
a/n  and  q/H.  For  practical  pixrposes,  the  dependence  of  Y  on  q/H  is  neg¬ 
ligible  when  q/H  is  sufficiently  larger  than  a/H.  The  relation  between  Y 
and  a/H  for  sufficiently  large  q/H  will  be  referred  to  as  the  K-callbration 
for  the  single-edge-notch  specimen  subject  to  splitting  forces.  The  K- 
calibration  resiilts  discussed  in  this  report  cover  a  range  of  values  of  a/H 
from  2  to  10,  which  is  sufficient  for  practical  purposes.  For  each  value  of 
a/H,  results  were  computed  for  values  of  the  parameter  q/H  of  a/H  +  1, 
a/H  +  2,  and  a/H  +  3  sufficient  to  define  the  practical  limit  of  influence  of 
this  parameter.  For  a/H  equal  to  4,  results  were  also  computed  for  q/H 
equal  to  a/H  +  6. 


AMALYTICAL  AMD  COMPUTATIONAL  PROCEDURE 

The  method  of  analysis  consists  in  finding  a  stress  function  X  that  sat¬ 
isfies  the  biharmonic  equation  V^X  =  0  and  also  the  boundary  conditions  at  a 
finite  niimber  of  stations  along  the  boundary  of  a  single -edge -notched  specimen, 
such  as  shown  in  figure  2.  The  biharmonic  equation  and  the  boundary  conditions 
along  the  crack  are  satisfied  by  the  Williams  stress  function  (refs.  8  and  9). 
Because  of  symmetry  (fig.  l),  the  coefficient  of  the  sine  terms  in  the  general 
stress  function  must  be  zero,  hence 

OO 

X(r,0)  =  ^  ^  |^_ 

n=l,2. . . 

,  f  T  -v  n .  n+1 
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The  stresses  in  terms  of  X  obtained  by  partial  differentiation  are  as  follows: 


,  S^x  sin  e  cos  0  4.  BX  sin^0 

(j  55  ss  ■  COS  V  -  Ci  '  ■■  ^  . 

y  ^  2  .  2  ^0  Sr  r  Sr  r 

ox  Sr 

^  SX  sin  9  cos  0  S^X  sin^0 

^  ^  S02 

a  =  sin^0  +  2  sin  9  cos  0  +  ^  cos^0 

^  Sy^  hr^  5^  Sr  ^  Sr  ^ 

2  ^  sin  0  cos  0  ^  S^X  cos^0 

S0  ^^2  ^2 

^  cos  20  S^X  sin  0  cos  9 

=  - -  sin  0  cos  0  — :r  +  -  r — rx  "  o  ^  9 

xy  Sy  Sr^  ^  37^  S0^ 

sin  0  cos  0  SX  _  cos  20  SX 
r  r^  B0 


(2) 


The  boundary-collocation  procediire  (ref.  5)  consists  in  solving  2m  si- 
miiltaneous  algebraic  equations  corresponding  to  the  known  values  of  X  and  ei¬ 
ther  ^X/5x  or  Sx/9y  at  m  selected  stations  along  the  boundary  ABCD  of 
figure  2j  thus,  values  for  the  first  2m  coefficients  of  the  Williams  stress 
function  are  obtained  when  the  remaining  terms  are  neglected.  Only  the  value 
of  the  first  coefficient  d^  is  used  for  the  present  purpose  since  the  stress - 

intensity  factor  is  equal  to  d^,  as  shown  in  reference  5. 

The  required  values  of  X  and  its  first  derivatives  at  the  m  selected 
boundary  stations  were  obtained  from  distributions  of  bending  moment  and  para¬ 
bolic  shear  stress  distribution  (fig.  2)  equivalent  to  the  actual  loads  shown 
in  figure  1.  The  equations  for  these  boundary  values  in  dimensionless  form 
are  as  follows : 


Along  A-B: 


Along  B-C: 


X  _  X  +  a  V 

P  B 


(3) 


Along  C-D: 
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The  distance  V  of  the  end  boundary  from  the  crack  tip  (fig.  2)  was  chosen 
to  be  0.8  a.  From  physical  considerations^  it  is  clear  that  the  boundary  should 
be  chosen  neither  close  to  the  crack  nor  close  to  the  load  point  because  of  the 
stress “field  disturbances  near  these  positions.  Preliminary  studies  similar  to 
those  described  in  reference  5  established  that  V/a  equal  to  0.8  was  about 
optimum^  and  that  minor  variations  from  0.8  had  a  negligible  effect  on  the  K- 
calibration- 

For  each  set  of  selected  values  of  the  primary  variable  a/H  and  the 
parameter  q/H^  the  collocation  computation  was  carried  successively  over  the 
range  of  m  equal  to  21  to  48  stations  in  increasing  increments  of  three.  Fig¬ 
ure  3  is  fairly  representative  of  the  resultant  variation  of  the  K-callbration 
against  the  number  of  boundary  stations.  Except  for  one  case^  the  variation  of 
the  dimensionless  stress -intensity  factor  Y  obtained  for  m  equal  to  42^  45^ 
and  48  was  less  than  1  percent  and  in  most  cases  it  was  a  small  fraction  of  a 
percent.  For  a/H  equal  to  2  and  q/a  equal  to  5/2^  the  variation  was  approx¬ 
imately  1^  percent.  From  this  observation  taken  together  with  the  nature  of  the 

trends  of  the  Y  values^  the  K-calibration  results  shown  in  table  I  were  based 
on  the  average  value  of  Y  for  m  equal  to  42^  45^  and  48. 


RESULTS  AND  DISCUSSION 

Table  I  shows  a  comparison  of  results  of  the  present  work  with  correspond¬ 
ing  results  from  reference  4  in  terms  of  the  dimensionless  quantity  KjBH^/^/ta 

as  a  function  of  relative  crack  length  a/H.  Over  the  range  of  relative  crack 
lengths  a/H  varying  from  2  to  10,  the  effect  of  relative  specimen  length  q/H 
becomes  negligible  when  not  less  than  a/H  +  2.  Figure  3  shows  a  typical  result 
in  evaluating  the  effects  of  a  large  q/H  ratio,  a/H  +  6,  on  the  dimensionless 
quantity  Y.  The  agreement  between  the  analytical  results  for  q/H  greater 
than  or  equal  to  a/H  +  2  and  the  experimental  results  of  reference  4  is  within 
1.3  percent  for  all  cases. 

For  a/H  less  than  2,  from  physical  considerations  one  can  conclude  that 
an  interaction  between  the  stress  field  of  the  applied  load  and  the  crack  tip  is 
likely  to  occur.  Thus,  use  of  simple  beam  theory  to  obtain  the  stress  function 
and  its  associated  derivative  at  location  V,  would  be  an  oversimplification  of 
a  complex  condition. 

The  plot  of  Y  against  the  reciprocal  of  relative  crack  length  (fig.  4)  is 
well  fitted  by  a  straight  line  passing  through  Y  =  3.46  at  H/a  =  0.  The 

value  3.46  is  equal  to  the  numerical  coefficient  obtained  as  a  first 

approximation  from  simple  beam  theory,  as  discussed  in  the  introduction.  Thus, 
extrapolation  of  the  collocation  results  to  the  limiting  case  of  a  very  long, 
slender  specimen  leads  to  a  value  which  is  in  good  agreement  with  the  elementary 
treatment. 
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MEASUREMENT  CAPACITY  IN  RELATION  TO  SPECIMEN  DEPTH 

The  ciirrently  accepted  criterion  for  evaluating  the  critical  value 

of  K  at  point  of  instability  of  crack  extension  in  first  or  open  mode^  is 
that  the  nominal  stress  at  the  crack  tip  should  not  exceed  the  yield  strength  of 
the  material  OyQ  in  a  valid  Kj^  test.  Thus^  the  maximum  nominal  stress 

6Pa/BH^  should  not  exceed  -which^  on  substituting  in  expression 

y2  _  gives  C^^y3(2H)  =  vhere  is  used  to  denote  the 

maximum  acceptable  value  Kj^.  Figure  5  shows  the  dependence  of  dimensionless 
Kj^  measurement  capacity  on  the  relative  crack  length  a/H  for  a  sufficiently 
large  effective  specimen  length.  The  specimen  measirrement  capacity  decreases 
with  increasing  relative  crack  length.  For  relative  crack  lengths  greater  than 
6^  the  drop  in  capacity  becomes  exceedingly  small. 

The  dependency  of  Kj^  measurement  capacity  on  the  specimen  thickness  is 
discussed  at  some  length  in  reference  1  and  is  unrelated  to  the  present  work. 

Lewis  Research  Center^ 

National  Aeronautics  and  Space  Administration^ 

Cleveland^  Ohio,  December  9,  1965. 
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APPEHDK  -  SYMBOLS 


a  effective  crack  length 


B  specimen  thickness 

C  estimate  of  maximum  value  of  that  can  be  measured  with  speci- 

IK  -Lc 

men  of  given  dimensions  and  yield  strength 


^n"'^2n-l 


coefficients  of  Williams  stress  function 


H 


K. 


K 


Ic 


m 

P 

q. 

r 

V 


Y 

G 


'X 


ys 


xy 


X 


depth  of  arm  of  split  beam 

stress-intensity  factor  of  elastic  stress  field  in  vicinity  of  crack 
tip 

critical  value  of  K  at  point  of  instability  of  crack  extension  in 
first  or  open  mode^  a  measure  of  plane  strain  crack  toughness  of 
material 

number  of  selected  boundary  stations  for  collocation  solution 

load  applied  to  both  ends  of  split  beam 

effective  specimen  length 

polar  coordinate  referred  to  crack  tip 

distance  from  crack  tip  to  boundary  for  collocation  analysis 
Cartesian  coordinates  referred  to  crack  tip 
dimensionless  stress -intensity  coefficient^  KBH^/^/pa 
polar  coordinate  referred  to  crack  tip 
stress  component  in  x-direction 

stress  component  in  y~direction 
0.2  percent  offset  tensile  yield  strength 
shearing  stress  component 
stress  function 
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TABLE  I.  -  COMPARISON  OF  RESULTS  OF  PRESENT  WORK 


WITH  CORRESPONDING  RESULTS  FROM  REFERENCE  4  IN 
TERMS  OF  THE  DIMENSIONLESS  QUANTITY  KjER^/'^/Pa, 
AS  FUNCTION  OF  RELATIVE  CRACK  LENGTH 


Relative 

crack 

length; 

a/H 

Results  of  " 

Boundary -collocation 
procedure 

Reference 

4 

( experimental 
results ) 

H  -  JJ  +  4 

i  -  ^  +  2 

H  H 

a  -  ^  +  3 

Y  «  K- 

j-BH^/yPa 

2 

5.88 

4.62 

4.65 

4.61 

3 

5.44 

4.25 

4.24 

4.21 

4 

5.22 

4.06 

4.05 

4.01 

5 

5.09 

3.94 

3.93 

3.90 

6 

5.00 

3.86 

3.85 

3.83 

7 

4.93 

3.81 

3.80 

3,78 

8 

4.90 

3.77 

3.76 

3.73 

9 

4.87 

3.73 

3.72 

3,70 

10 

4.82 

3.71 

3.72 

3,67 

Figure  1.  -  Edge- cracked  plate  specimen  subject  to  splitting  forces. 


Number  of  boundary-collocation  stations,  m 

Figure  3.  -  Value  of  dimensionless  stress-intensity  coeffi¬ 
cient  Y  against  number  of  boundary-collocation  stations 
for  relative  crack  length  a/H  =  4. 


Dimensionless  stress-intensity  coefficient, 


Figure  4.  -  Dimensionless  stress-intensity  coefficient  as 
function  of  reciprocal  of  relative  crack  length. 


Figures.  -  Dependence  of  measurement 
capacity  of  single-edge-notched  speci¬ 
mens  subjected  to  splitting  forces  on  relative 
crack  length. 


NASA- Langley,  1966  E-3271 


aeronautical  and  space  activities  of  the  United  States  shall  he 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof.” 

— National  Aeronautics  and  Space  Act  of  1958 
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